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Using state-of-the-art first-principles calculations we study the magnetic behaviour of CeOFeAs. 
We find the Ce layer moments oriented perpendicular to those of the Fe layers. An analysis of 
incommensurate magnetic structures reveals that the Ce-Ce magnetic coupling is rather weak with, 
however, a strong Fe-Ce coupling. Comparison of the origin of the tetragonal to orthorhombic 
structural distortion in CeOFeAs and LaOFeAs show marked differences; in CeOFeAs the distortion 
is stabilized by a lowering of spectral weight at the Fermi level, while in LaOFeAs by a reduction 
in magnetic frustration. Finally, we investigate the impact of electron doping upon CeOFeAs and 
show that while the ground state Fe moment remains largely unchanged by doping, the stability 
of magnetic order goes to zero at a doping that corresponds well to the vanishing of the Neel 
temperature. 

PACS numbers: 74.25. Jb,67.30.hj, 75.30.Fv,75.25.tz, 74.25. Kc 



The recently discovered [J, l^l family of FeAs based 
compounds that, upon electron doping, become super- 
conducting with transition temperatures up to 55K, are 
attracting a lot of interest. Structurally these materi- 
als ROFeAs (R=La, Ce, Pr, Nd, Sm) are very similar 
in that they are formed from FeAs layers separated by 
rare earth or lanthanide oxide layers. In striking contrast 
to the well known cuprates, these materials are weakly 
magnetic metals, leading to the possibility of large spin 
fluctuation effects In particular, at the onset of 

superconductivity the moment on the Fe atoms vanishes, 
and a key question concerns the possible role of such spin 
fluctuations in the superconducting transition]^, 01 • 

Despite the diverse set of rare earth and lanthanide 
atoms involved in these materials, physically they share 
many similarities. In particular, (i) at temperatures 
around ~150K a structural phase transition from tetrag- 
onal to orthorhombic crystal symmetry occurs and, (ii) 
this is then closely followed by a magnetic phase tran- 
sition to a spin order anti-fcrromagnetic (AFM) in na- 
ture, (iii) upon doping with Florine the AFM order is 
suppressed and superconductivity appears. Furthermore, 
calculations reveal that (iv) the non-magnetic Fermi sur- 
faces are all strongly nested[l, 0, Q, and (v) that the 
moment of the Fe atoms depends critically upon the sep- 
aration of the Fe layer from the adjacent As layer[l,0,S]. 

Given the diversity of constituents involved in this class 
of materials, such uniformity of behaviour appears, at 
first sight, somewhat surprising. The question then arises 
if the underlying physical mechanisms behind such phe- 
nomenon as the structural distortion or doping behaviour 
are also the same. In the present work we investigate 
this by analyzing the properties of CeOFeAs as com- 
pared to the well studied LaOFeAs. The choice of these 
two materials is motivated by the fact that amongst the 
rare earth and lanthanide oxypnictides they show rela- 



tively large differences for two key physical properties; 
(i) the low temperature moment of the Fe atom, 0.35/iB 
in LaOFeAs is the lowest while in CeOFeAs (0.94ub) 
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is the highest amongst all oxypnictides [9|, ll( 
and (ii) the superconducting transition temperature also 
differs strongly; 26 K in LaOFeAs [E ^ and 45 K in 
CeOFeAs [3 [lllii. 



Remarkably, we find that the mechanism behind the 
structural phase transition in CeOFeAs and LaOFeAs 
is quite different; the former case being driven, essen- 
tially, by the one electron kinetic energy i.e. Fermi 
surface effects, while the latter is driven by magnetic 
frustration^, [l3|. In addition, the behaviour upon elec- 
tron doping is substantially different. We find that for 
CeO (i_a;)Fa;FeAs the moment is almost unchanged upon 
electron doping (a reduction of 3% for x — 0.10), whereas 
calculations of LaO(i_2,)Fj;FeAs have revealed a strong 
suppression of the moment (90% at a; = 0.10). This may 
be reconciled with the vanishing of the Neel temperature 
at x = 0.06 by the fact that the magnetic order becomes 
meta-stahle near this doping. 

In the present work all calculations are performed using 
the state-of-the-art full-potential linearized augmented 
plane wave (FPLAPW) method fl8^, implemented within 
the Elk code 19]. To obtain the Pauli spinor states, the 
Hamiltonian containing only the scalar potential is diag- 
onalized in the LAPW basis: this is the first-variational 
step. The scalar states thus obtained are then used 
as a basis to set up a second-variational Hamiltonian 
with spinor degrees of freedom ^18'|. This is more effi- 
cient than simply using spinor LAPW functions, but care 
must be taken to ensure that there is a sufficient num- 
ber of first-variational eigenstates for convergence of the 
second-variational problem. We use a shifted k mesh of 
10 X 10 X 6, and 260 states per k-point which ensures 
convergence of the second variational step as well as the 
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convergence with respect to the k-points. All the exper- 
imental lattice parameters from Ref. Q have been used. 



as indicated in experiments 1J| we find that the Ce and 
Fe moments are strongly coupled; a small change (1%) of 
the Ce moment (performed with a fixed spin calculation), 
leads to a large change (3%) in the Fe moment. 
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FIG. 1: (color online) Top-left panel: energy (in meV) per 
formula unit as a function of position of the As atom cal- 
culated using LDA and GGA. Bottom- left panel: same as 
upper panel but the calculation is spin-polarized. Right-hand 
panel shows the moment per Fe atom (in fis) calculated using 
LSDA and GGA. 
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FIG. 2: (color online) Top panel shows the energy (in meV 
per formula unit) and the middle and bottom panels show the 
magnetic moment (in /is) per Ce and Fe atom respectively. 
All quantities are plotted as a function of the spin-spiral q- 
vector. The inset is a schematic of the spin configuration of 
the Fe and Ce atoms. 



We first consider the undistorted undoped ground state 
of CeOFeAs. One of the most striking aspects of the 
theoretically intensively studied compound LaOFeAs has 
been the spread of results for magnetic properties, at- 
tributed to both an unusual sensitivity to the approxi- 
mation to exchange-correlation, and a sensitive depen- 
dence upon the separation of the Fe and As layers, 
zas 0, 0,S| ■ This latter behaviour is also found in the case 
of CeOFeAs [20|. however in contrast to LaOFeAs a spin 
polarized generalized gradient approximation (GGA) [2l| 
calculation yields an optimized zas in near perfect agree- 
ment with experiment, see Fig.[T] 

Given this choice of zas we now determine the 
ground state magnetic structure using, in addition to the 
GGA functional, the loca l sp in density approximation 
(LSDA)!!!] and LSDA-LUjli functionals (for the latter 
we use U = 6 eV and J = 1 eV for the Ce atom). Con- 
sidering first only coUinear structures, we find that the 
Fe layer adopts the stripe AFM structure found in all the 
oxypnictides, with the magnetic moments of the Ce layer 
perpendicular to those of the Fe layer. The Fe (Ce) mo- 
ments are found tobeI.58/i_B (0.56/x_b), {0.92^b), 
and l.SOfiB (0.54/is) for the LSDA, LSDA-LU, and GGA 
functionals respectively, with the corresponding experi- 
mental values 0.94/^B for Fe and 0.83/xb for Ce. The 
agreement with experiment for the Fe moment is thus 
rather poor for all functionals considered. This may be 
attributed to the exceptionally sensitive dependence of 
the Fe moment on zas, see Fig. [1] which for CeOFeAs is 
even more pronounced than in LaOFeAs. Interestingly, 



In order to study the possible incommensurate spin 
structure we have calculated the total energy as a func- 
tion of the spin-spiral vector q, for various directions in 
the Brillouin zone (BZ). For the in plane spin spirals 
(going in the direction [0.5, 0.5, 0] to [1, 1, 0]) we find 
a clear sharp minimum at the commensurate q vector of 
[0.5, 0.5, 0], equivalent to stripe AFM spin configuration 
(top panel Fig. [J). Interestingly, however, we find that 
the various spin configurations are almost degenerate in 
the direction [0.5, 0.5, 0.4] to [0.5, 0.5, 0.5]. This is in- 
dicative of weak inter-plane Ce-Ce and Fe-Fe coupling; 
this findin g is concomitant with experimental results of 
Zhao et al. 
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Turning now to the spin spiral moments (lower two 
panels Fig. [J), we find that the moment remains almost 
unchanged upon changing the magnetic structure away 
from the stripe AFM spin configuration. This is in strik- 
in g co ntrast to LaOFeAs, which remains spin polarized 
[sl l24i[25t only in a small region about the M point ([0.5, 
0.5, 0]), with the moment vanishing elsewhere in the BZ. 
Since for CeOFeAs the moment on the Fe atoms changes 
only slightly - 5.2% - upon moving across the BZ hence 
the role of spin fluctuations could be rather different in 
these two materials. 

This difference may be understood as a consequence 
of the nature of the various mter-plane and mira-plane 
magnetic couplings in this material. In particular, the 
infra-plane coupling of the Ce atoms is very weak - we 
find the energy difference between FM and AFM ordered 
in-plane Ce moments to be almost degenerate with, ad- 
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ditionally, the magnitude of the Ce moments unchanged 
by this choice of FM or AFM order. Thus the nature 
of the Ce-Ce interaction results in the Ce moment re- 
maining unchanged by a spin wave configuration; this in 
turn acts to preserve the Fe moment due to the strong 
inter-plane Fe-Ce coupling. 



to the frustration driven mechanism of LaOFeAs. In this 
regard in Fig. [4] are shown the non-magnetic as well as 
magnetic (AFM stripe phase) density of states (DOS) for 
the undistorted and distorted lattices. One feature coni- 
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FIG. 3: (color online) Top-left panel shows the total en- 
ergy and bottom-left panel the one-electron kinetic energy (in 
meV) per formula unit. The right panel shows the moment 
per Ce (top) and Fe (bottom) atom (in fis) as a function of 
distortion angle (in degrees). 

As with all the oxypnictides, CeOFeAs undergoes 
a structural phase transition from tetragonal to or- 
thorhombic crystal symmetry IJ] , which in this com- 
pound occurs at 160K. In order to determine the physical 
reason behind this transition we have performed ab-intio 
LSDA calculations. As may be seen in Fig. [3] the non- 
magnetic compound does not show any crystal distortion, 
but upon performing spin polarized calculations one finds 
a minimum at a distortion angle of 90.30°. The differ- 
ence in energy between the non-magnetic undistorted and 
magnetic distorted materials is 17meV (197 K), which is 
comparable to the experimental critical temperature of 
160K. Thus both the distortion and critical temperatures 
are very similar to those found in LaOFeAs however, as 
we will now show, the underlying mechanism of the tran- 
sition is entirely different. 

In the case of LaOFeAs the crystal distortion could be 
rationalized as being driven by an increase in moment 
upon distortion [sj, due to a lowering of magnetic frustra- 
tion in the distorted lattice [1, [l3]- On the other hand, 
the moment of both Fe and Ce atoms remains essen- 
tially unchanged in CeOFeAs upon distortion, see right 
hand panels of Fig. [31 thus ruling out such a mechanism 
for this compound. However, the one-electron kinetic en- 
ergy of CeOFeAs shows a clear minimum at 90.30° while, 
in dramatic contrast, the LaOFeAs one electron kinetic 
energy monotonically increases as a function of distor- 
tion angle. Thus the structural distortion in CeOFeAs is 
driven by an electronic structure mechanism, in contrast 
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FIG. 4: (color online) Density of states calculated using 
the LSDA for undistorted non-magnetic, undistorted mag- 
netic and distorted magnetic structures for CeOFeAs and 
LaOFeAs. Stripe anti-ferromagnetic order was used for the 
magnetic calculations. 

mon to both compounds is that, as one would expect, 
spin polarization results in a large shift of spectral weight 
away from the Fermi level. Turning to the impact of dis- 
tortion one observes a striking difference: in LaOFeAs 
there is a substantial increase in spectral weight at the 
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26l | while, in contrast, in CeOFeAs spec- 
tral weight is moved away from the Fermi level. Thus 
in this latter case the redistribution of spectral weight 
drives the transition, while for LaOFeAs the opposite is 
the case. Clearly, therefore, although both the distor- 
tion angle and transition temperatures are very similar 
in these two oxypnictides, the origin of the structural 
phase transitions is very different. At this point it is 
also worth mentioning that heat capacity measurements 
of McGuire et al. showed that this decrease of the DOS 
at the Fermi level in CeOFeAs on distortion is consistent 
with the Seebeck coefficient measurements [27j . 

The most important property of the iron oxypnictides 
is the occurrence of a superconducting phase transition 

An inter- 



at a critical electron doping[lJ, [1 



esting difference between the temperature-doping phase 
diagrams of CeOi_a;Fa;FeAs and LaOi_a;Fj;FeAs is that, 
for the former case, the Neel temperature of the mag- 
netic phase goes continuously to zero [3] as critical dop- 
ing is approached {xc = 0.06), while in LaOi_a;Fa:FeAs 
one instead finds a sharp dropjl, [sli in the Neel tem- 
perature at a critical doping of Xc — 0.045 (see right 
hand panel Fig. [5]). Concomitantly, low temperature 
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FIG. 5: (color online) Left panel: moment per Fe atom (in /is) 
as a function of doping calculated using the LSDA, LSDA+U 
and GGA. The right panel: shows the difference between the 
magnetic and the non-magnetic total energy (in K) calculated 
using the LSDA. The experimental data for CeOi-xF^FeAs 
are taken from Ref. Q and for LaOi_a;Fa;FeAs from Ref. 



measurements of the magnetic order of the Fe layer in 
CeOi_2;Fa;FeAs show that it is entirely lost before Xc (left 
panel Fig. [5]). An important question for understanding 
the role of magnetism in the superconducting transition 
is then whether it is the moment of the Fe atoms that 
vanishes, or whether it is simply the stripe AFM order 
that vanishes. 

To investigate this issue we calculate CeOi_a;F^FeAs 
by deploying the virtual crystal approximation (VGA). 
Although the VGA neglects much of the physics of dis- 
order, it has been shown recentlvfs^ to provide a sur- 
prisingly accurate account of magnetism and Fermiology 
in La0i_2;F2:FeAs. In order to determine the impact 
of the approximation to the xc functional, calculations 
were performed using the LSDA, LSDA+U, and GGA 
functionals. 

Amongst the undoped Fe oxypnictides the Fe moment 
is largest in GeOFeAs and smallest in LaOFeAs, and it 
is therefore surprising that in the latter compound the 
magnetic phase appears more stable for x < Xc- In- 
deed, calculations of the Fe moment in GeOi-^F^FeAs 
(left panel Fig. [5]) show that, in apparent contradiction to 
the experimental data, there is almost no change upon 
doping. However, if we consider the magnetization en- 
ergy, i.e. the quantity i^non-mag - £^mag (right panel 
Fig. [5|), this falls sharply at exactly the critical doping 
Xc = 0.06. Thus it is the stability of the magnetic order 
which is reduced upon doping, not the actual Fe moment. 
This sudden reduction of stability, which is not found in 
LaOi_a;F3;FeAs, is behind the differing forms of the Neel 
temperature phase boundary. 

In the case of LaOFeAs calculations have shown 8] 
that, upon doping, the distortion angle changes from an 
undoped value of 6^ = 90.27° to 6* = 89.85°; addition- 



ally, there is the appearance of incommensurate mag- 
netic phases. We have investigated if such phases appear 
in the doping phase diagram of GeOi_a;Fa;FeAs and find 
that (i) the distortion angle is unchanged upon doping 
(0 < a; < 0.06) and (ii) the ground state magnetic struc- 
ture remains stripe AFM. 

To conclude, we have performed a comparative 
study of the two oxypnictides GeOi_a;Fa;FeAs and 
LaOi_a;Fa;FeAs. Although many common features, e.g. 
a similar non-magnetic Fermi surface, may be found in 
the iron oxypnictides we have shown here that profound 
differences also exist. In particular we find that the struc- 
tural distortion in GeOFeAs is driven by one electron ki- 
netic energy (i.e. electronic structure effect), in contrast 
to LaOFeAs where a lowering of magnetic frustration has 
been identified as the mechanism. Furthermore, the be- 
haviour under electron doping differs markedly between 
the two materials; in GeOi_a;Fj;FeAs we find that while 
the ground state moment is essentially unchanged with 
doping, the stability of the moment is sharply reduced. 
Finally, we have - via calculations of incommensurate 
spin spiral structures - carefully investigated the vari- 
ous magnetic couplings of the constituents of GeOFeAs; 
we find that the Ge atoms are only weakly coupled to 
each other, while the Fe atoms are much more strongly 
coupled to each other as well as to the Ge. 
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